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Abstract. Action at distance in Newtonian physics is replaced by 
finite propagation speeds in classical post-Newtonian physics. As 
a result, the differential equations of motion in Newtonian physics 
are replaced by functional differential equations, where the delay 
associated with the finite propagation speed is taken into account. 
Newtonian equations of motion, with post-Newtonian corrections, 
are often used to approximate the functional differential equations. 
Are the finite propagation speeds the origin of the quantum me- 
chanics? In this work a simple atomic model based on a functional 
differential equation which reproduces the quantized Bohr atomic 
model is presented. As straightforward application of the result 
the fine structure of the hydrogen atom is approached. 



1. Introduction 

Newtonian forces (for example, the inverse square law for gravita- 
tion) imply "action at distance". This absurd, but outstandingly suc- 
cessful, premise of Newtonian theory predicts that signals propagate 
instantaneously. In classical physics, relativity theory postulates that 
signals propagate with a velocity that does not exceed the velocity of 
light. Thus, the forces of Newtonian physics must be replaced by force 
laws that take into account the finite propagation speed of the classical 
fields which determine the forces acting on a moving body. In turn, the 
ordinary or partial differential equations of Newtonian physics, which 
are derived from the second law of motion mr = F, must be replaced 
by corresponding functional differential equations where the force F is 
no longer a function of just position, time, and velocity; rather, the 
classical force law must take into account the time delays due to the 
finite propagation speed of the classical fields. 
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The functional differential equations of motion for classical field the- 
ory are generally difficult, often impossible, to express in a form that is 
amenable to analysis. Thus, in order to obtain useful dynamical pre- 
dictions from realistic models, it is frequently to replace the functional 
differential equations of motion by approximations that are ordinary 
or partial differential equations, see The purpose in these works is 
to discuss some of the mathematical issues that must be addressed to 
obtain a rigorous foundation for the post-Newtonian dynamics, that is, 
Newtonian dynamics with relativistic corrections, see for instance |3] 
and the references therein. For the electromagnetic classical field, in the 
ideal case of a point-charge particle, the resulting retarded potentials 
are the Lienard-Wiechert potentials. For the gravitational classical 
field we must use the Einstein's field equation. Simple models of these 
equations are the subject of current research. The basic idea of post- 
Newtonian approximation, from a mathematical point of view, is the 
expansion of model equations in powers of 1/c. From a physical point 
of view, the idea is to consider low velocity (compared with the speed of 
light). Note, for example, that the relativistic form of Newton's second 
law, where the rate of change of the momentum is given by 



reverts to Newton's law in the low-velocity limit. 

According to Maxwell's field equations, a charged particle produces 
electromagnetic fields as it moves. Since, in this particle radiates 

energy, it must slow down. In the theory of the electron point charge 
by considering motion along in a line, Dirac propose a self-force (the 
radiation reaction force) given by 



which is the half difference of the retarded and advanced forces, where 
q is the charge of the electron, see [TT]. Therefore, a post-Newtonian 
model for the motion of an electron, confined to move on a line and 
with radiation reaction taken into account, is given by the Abraham- 
Lorentz equation 



where F is an external force. Since the electron radiates (produces fields 
that carry energy) the self force should cause the particle to lose energy 
and slow down. That's why, the presence of the third derivative term 
in the first differential equation is called radiation damping. However, 




(1) 



mi = — -'x + F, 



c 
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in these post-Newtonian models (where the differential equations are 
not of second order) the "runaway" solutions appear, see ^T]. For 
instance, in absence of external forces, equation (JTJ) reduces to 

2q 2 ... 

mx = — - x , 
3c d 

and this equation has the solution x — C where C is an arbitrary 
constant, and other solutions where the acceleration is proportional 
to exp(3mc 3 t/(2g 2 )). Hence, the acceleration grows indefinitely with 
time. This means that, a charge which goes out of a field, when leaving 
it, must self- accelerate indefinitely; which is an absurd. These runaway 
solutions are clearly not physical. What do they represent? How should 
they be eliminated? and What is the correct Newtonian equation with 
the radiation damping taken into account? 

The mathematical answer to all these questions is approached in J3] 
and the subsequent works where these post-Newtonian models 

are recognized as singularly perturbed Newtonian equations. In order 
to recover the correct Newtonian equations with the post-Newtonian 
corrections, the Fenichel's geometric singular perturbation theory is ap- 
plied (in particular, the reduction to the slow-manifold). These New- 
tonian equations with the post-Newtonian corrections give physically 
reasonable dynamics; in particular, the runaway solutions are elimi- 
nated. Anyway, how can we justify using these models? Note, for 
instance, that the slow-manifolds in these models are unstable; nearby 
runaway solutions. In applied mathematics, we usually justify approxi- 
mations by their stability. To validate the slow-manifolds reductions it 
must be shown that the resulting Newtonian model equations are sta- 
ble with respect to the dynamics of the original functional differential 
equations, the true equations of motion in classical physics. Therefore 
further investigations are required in this direction for the study of the 
delay equations. 

However, it is interesting to note that the presence of a small delay 
in a conservative system often results in damped long-term dynamics 
on an associated inertial manifold, see 011]. For example, the Duffing- 
type model equation 

x + u 2 x = —ax(t — t) + bx 3 (t — t), 

with small delay r in the restoring force, reduce (by a formal compu- 
tation to first order in r) to the van der Pol-type model equation 

x, + r(3bx 2 — a)x + (a + uj 2 )x — bx 3 = 0, 

on its inertial manifold. As it has been noticed in -3] this example 
illustrates a phenomenon that is a reminiscent of quantization: while 
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most periodic solutions in one parameter families of periodic solutions 
in a conservative system disappear in the presence of a small delay, 
some persist as limit cycles. The author of jSj asks himself whether 
this observation has a physical significance. 

The solutions of the functional differential equations can, however, 
admit an infinite discrete spectrum. For example, we consider the 
retarded harmonic oscillator, given by the linear, second order, retarded 
functional differential equation 

x + x(t — r) = 0, 

with small delay r. In ^H] it is showed that this equation exhibits 
an infinite spectrum of discrete frequencies. Its general solution is a 
convergent linear combination of oscillations at an infinity of discrete 
("quantized") frequencies. As in quantum mechanics, in order to de- 
termine a unique solution an initial function needs to be known. The 
above consideration remains valid for any linear functional differential 
equation with constant coefficients and constant delay. Moreover, the 
locally linear approximation suggests that such "quantization" is also 
to be expected for non-linear functional differential equations. 

At the beginning of the 20th century, Planck [14] initiated the quan- 
tum mechanics with his contribution to the black body radiation. Ein- 
stein jH], following the ideas of Planck [T^j, contributed to the devel- 
opment of the theory of quanta which is the embryonic step needed to 
arrive to the Quantum physics theory. It is interesting to note that 
Poincare had been implicated in the discussion of the quantum theory, 
but the premature death of Poincare deprives of its contributions in 
this theory, see for instance fUJEl- in fact, Poincare participated in 
the first congress of Solvay in October of 1911 and he died in July of 
1912. One of the interesting known phenomenon studied by Poincare is 
the concept of limit cycle, see [T^]. This phenomenon does not occur in 
the Hamiltonian systems studied by standard physics theories. It only 
appears in systems with friction, i.e., systems with dissipative energy. 

It is also interesting to note that one of the problems which orig- 
inates the quantum theory was the problem that appears when the 
idea of the planetarium system is applied to the atomic model. This 
idea was proposed by Rutherford in 1911 20 as a consequence of the 
experimental results obtained when bombing an atom with a parti- 
cles. The problem in the model of Rutherford was that the charged 
electrons are accelerated in their movement around the nucleus and by 
the electromagnetic classical theory any accelerated charged body ra- 
diates energy. The problem of the atomic stability was initially solved 
by Bohr in 1913 [1 and it marks the success of the quantum theory 
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and its posterior development. The atomic model of Bohr predicts the 
radiation spectrum of certain atoms and the quantization of the energy 
in different energy levels is then obtained. 

If you see the development of the quantum theory from the initial 
contributions, it is evident that each step is made with extra assump- 
tions. For instance, the introduction of the quanta in the radiation of 
a black body and in the foto-electric effect, cf. [S], the quantization of 
the energy in the movement of an electron which moves as an harmonic 
oscillator under the influence of an harmonic restoring force, cf. |14j . 
Another example is the quantization of the angular orbital impulse of 
an electron in an atom, although the electron in an atom is accelerated 
in its movement around the nucleus. In this last case, it is assumed 
that this electron does not radiate energy, see pQ. However, we notice 
the difference between the Bohr quantization of the angular orbital im- 
pulse of an electron, which moves under the Coulomb force (L = nh 
for n = 1, 2, 3, . . ., where h is a multiple of the Planck constant h di- 
vided by 27r), and the Planck quantization of the energy of a particle, 
as an electron, which moves as an harmonic oscillator (E = nhv for 
n — 1, 2, 3, . . .). In fact the quantization of the angular orbital impulse 
of an electron leads to the quantization of the total energy but with an 
equation quite different than the Planck equation. 

In jH] it is showed that the intrinsic phenomenon (the quantization of 
the energy) that appears in the first and simple systems initially studied 
by the quantum theory as the harmonic oscillator and the movement 
of a charged particle under the Coulomb force, can be obtained from 
the study of dissipative systems. In other words, it is showed that 
the phenomenon of the quantization of the energy of a particle which 
moves as an harmonic oscillator can be obtained via a classical system 
of equations. The same assertion also applies to the phenomenon of 
the quantization of the energy of a charged particle which moves under 
the Coulomb force and which loses and wins energy (for example if 
we consider the classical case where the electron radiates and absorbs 
energy from the electric field of the nucleus). Therefore, these phenom- 
ena are not intrinsic of the quantum theory, but also appear in classical 
systems. In fact, they appear in the qualitative theory of differential 
equations developed by Poincare from 1881 [To] . 

Nevertheless, the most important problem is to find the exact form 
of the dissipative term and the interpretation of its physical mean- 
ing, see jH]. The retarded case, already explicitly incorporates certain 
subtle mathematical features of electrodynamics and relativity noticed 
by Poincare, but overlooked by Einstein and subsequent researchers. 
Based on the study of the retarded systems, a simple atomic model 
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given by a functional differential equation which reproduces the quan- 
tized Bohr atomic model is presented in this paper. 

The paper is not at all an alternative to the quantum theory, because 
the large development of the quantum theory in all the past century, 
the success in all its predictions, is outside of all doubt. This work 
does not pretend, in any case, to substitute quantum mechanics but 
to complete the knowledge that it gives to us. On the other hand, the 
proposed model does not reflect the whole rich behavior of the quantum 
modern theories developed from 1925 by Schrodinger [2H 1221 123] > Born 
[2], Heisenberg [TU], Dirac jjj, and others. The goal of the paper is 
to ask if the finite propagation speeds is the origin of the quantum 
mechanics. 

To begin with, in [18j . it is assumed that the two particles are ro- 
tating rigidly in circular orbits around a common center of masses. 
Moreover, a force which varies inversely as the square of the retarded 
distance is considered. The retarded distance is the distance from the 
current position of the electron to the "last seen" position of the proton. 
The simple expression for the force helps us to intuitively understand 
the consequences of a delay and under such circumstances, the angular 
momentum cannot be conserved. Thus, we have the astonishing situa- 
tion that, purely under the action of internal forces, the system suffers 
a net torque. Now, the radiation term is introduced. But, the exact 
form of the radiation damping term is not clear. Finally, the simple 
heuristic case of the retarded inverse square force is used, to determine 
whether there can be a balance of forces between the delay torque and 
the 3rd order radiation damping. A total success does not take place 
and a value of r which is smaller than the Bohr radius is obtained. 
Nevertheless, the author, in [TBI , affirms that further investigations are 
required to determine the exact effects of radiative damping, and that 
it was prematurely concluded that radiative damping makes the clas- 
sical hydrogen atom unstable. In the following section we present a 
simple atomic model based on a functional differential equation which 
reproduces the quantized Bohr atomic model. It is important to stand 
out that what we will carry out in the following section is not a post- 
Newtonian approach in which r is small. This is what has been made 
up to now and in the mentioned works El El El- From now on, we 
will accept that the laws governing the movement have a delay (a de- 
lay that does not need to be small) and we will find a solution of the 
functional differential equation in a very simple case. 
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Figure 1 . The retarded electrodynamic 2-body problem. 

2. The retarded electrodynamic 2-body problem 

We consider two particles interacting through the retarded inverse 
square force. The force on the electron exerted by the proton is given 
by 

(2) F = r. 

The force acts in the direction of the 3-vector r, along which the proton 
is " last seen" by the electron. The 3-vector r may be represented by 

r = r p (t - r) - r e (t), 

where r p (t) and r e (t) denote respectively the instantaneous position 
vectors of the proton and electron, respectively, at time t, and r is 
the delay, so that r p (t — r) is the "last seen" position of the proton. 
Assuming that the two particles are in rigid rotation with constant 
angular velocity cu, and referring back to Fig. 1, we have, in 3-vector 
notation, 

r e = r\ [cos out i + sin uit j] , 

and 

r p = — r 2 [cos uj{t — r) l + sincj(t — r) j]. 
Hence, the 3-vector r is given by 

r = [— r 2 cosuj{t — t) — r\ coswt] l + [— r 2 sinw(t — r) — r\ sin cut] j, 

Now, we introduce the polar coordinates (r, 6) and define the unitary 
vectors 1 = cos 9 i+sin 9 j and n = — sin 9 i+cos 9 j. By straightforward 
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calculations it is easy to see that the components of the force (J2J) in the 
polar coordinates are 



e 2 e 2 



F r = K—r ■ 1 = (— r 2 cos(W) - ri)K— 

and 

e 2 e 2 
F e = K—r ■ n = r 2 sinfcj-r).^— 

The equations of the movement are 

(3) mf — mrO 2 = F 7 

(4) mrO + 2mr6 = F e . 

The second equation (0J) can be written in the form 

(5) --77 = -^r{rnr 2 e) = F e = r 2 sin(W)K^r. 

r at r at r 6 

In 1913 Bohr pQ introduced the quantization of the angular momen- 
tum of the form L = nh/(27r) where h is the Planck constant. If we 
accurately study equation (0) we see that the analytic function sin(u;r) 
has a numerable number of zeros given by 

(6) UJT = kn , 

with fceZ, which are stationary orbits of the system of equations Q 
and (0J). When out ^ kn we have a torque which conduces the electron 
to the stationary orbits without torque, that is, with out = kn. 

This is a new form of treating the hydrogen atom from a dynamic 
point of view instead of from a static point of view, as it has been made 
up to now. Moreover, in this model the delay r is not small, in fact 

kir kn kT 



u 2tt/T 2 

where T is the time taken by an electron to complete its orbit. There- 
fore, the delay is a multiple of the half-period T/2. 

On the other hand, in a first approximation, the delay r can be 
equal to r/c (the time that the field uses to goes from the proton to 
the electron at the speed of the light). In this case, from equation © 
we have 

kn r 

(7) t = — = - . 

LO c 

Taking into account that to = vg/r, from (J7|) we have vg/c = kir. 
However, from the Relativity theory we know that vg/c < 1, then we 
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must introduce a new constant g in the delay. Hence, r = gr/c and 
the new equation (J7J) is 

(8) T= k 1=9 _r^ 

uj c 

and now Vg/c = kn/g, i.e. vg = kirc/g and from (jHJ) we also have 
r = kirc/(guj). In our model case of a classical rigid rotation we have 
9 = ujt with uj > 0. Therefore, 9 = uj and 9 = 0. Hence, equation Q 
for ujt = kir is 

2mruj = 0, 

which implies f = and r = where is a constant for each A;. On 
the other hand, equation (j2J) for ujt = kir takes the form: 

,2 „2 J2 



(9) -mr# 2 = -m^ = (-r 2 (-l) n - ri)i^ w -rK-^, 

assuming that r ~ r x due to r 2 C rj. 

From the definition of angular momentum L = mr 2 9 = mr 2 uj = 
mrvg we have that v$ = L/(mr). Substituting this value of Vg into 
equation Q we obtain r = L 2 /{mKe 2 ). The energy of the electron 
(substituting the values of vg and r) is given by 



mv 



2 „2 



(10) E = — °--K 



e 2 K 2 m e 4 



2 r 2L 2 

The angular momentum for ujt = kir is 

kirc L 2 

(11) L = mvgr = m- 



g mKe 2 ' 

which is an equation for the angular momentum. Isolating the value 
of L we obtain L = (Ke 2 g)/ (knc). If we introduce this value of the 
angular momentum in the expression of the energy (|1U|) we have 

2 2 

^ mix c , 9 

E = —k 2 , 

2g 2 ' 

In 1890 Johannes Robert Rydberg generalized Balmer's formula and 
showed that it had a wider applicability. He wrote his formula as 

A R \n\ 7t|) 

where A is the wavelength, ni and ri2 are integer numbers and R is 
known as the Rydberg constant. Therefore, the energy levels are pro- 
portional to 1/n 2 and, of course, negative, because these are bound 
states, and we count energy zero from where the two particles are in- 
finitely far apart. Hence, (identifying n = \k\) we must impose that 
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the constant g = k 2 gi where g\ is another new constant and then the 
energy takes the form 

rmr 2 c 2 1 

E = 



2g 2 k 2 ' 

If we recall the expression of the energy levels given by Bohr in 1913 

4 

E 



%e 2 h 2 n 2 ' 

we can compare these two expressions of the energy levels and we get 
the explicit value of the Planck constant 

(12) h eV 



2tc e c 

The dimensional analysis gives [J ■ s] = [C 2 ]/([C 2 /(Jm)][m/s]) which 
is correct. The introduction of a new fundamental constant is avoided 
(as it happens in the quantum mechanics with the Planck constant) 
because through the delay the speed of the light c appears. The ap- 
pearance of this dimensional constant c, usually absent in the non- 
relativistic quantum mechanics, allows to give the expressions of the 
physical quantities with the correct dimensions. Substituting into equa- 
tion (|T2*Jl the value of the electron charge e, the electric permittivity 
constant Eq and the value of the speed of light c we obtain that the 
adimensional constant gi must take the value g\ = 429.868. 

We notice that this value of g\ = 429.868 is the inverse value of 
the fine structure constant \ja = 137.036 multiplied by ir, that is, 
<7i = ir /a. Moreover, the value of this fine structure constant a is 

- 

and when substituting this value of gi in the explicit expression of the 

Planck constant (|T2"|) the equality is identically satisfied. Therefore we 

have found the value of the adimensional constant g\ and consequently 

the expression of the delay r which is 

or ,t)Qir , 9 iir 
13 t = — = k — = k — . 

c c ac 

From the found value of the angular momentum and the value of 
vq = ac/k we have 

, . _ Ke 2 g Ke 2 k ac 

(14) L = — = = mvgr = m—r. 

kite ac k 

Isolating the value of r from equation ()14j) we obtain 

Ke 2 k 2 
r = . 

ma 2 c 2 
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Taking into account that K = 1/ (Atteo) and the value of the fine struc- 
ture constant a, we arrive to the classical radii of the stationary orbits 



(15) r. = ^4. 

Moreover, the relation vg / c = kit / g = a j k is consequent with the defi- 
nition of a. One of the interpretations of the fine structure constant a 
is that a relates the speed of the electron in the lowest energy level in 
the atom of hydrogen with the speed of the light. This is straightfor- 
ward because if we substitute the expression (JT5|l of r n in the classical 
expression of Vg = nh/(2irmr n ) given by the quantum mechanics, for 
the case n — 1, and we divide by the speed of the light c, we obtain 
that vg/c = a. 

Summarizing, we could have begun our analysis with the found delay 
definition (|13j) . because our model reproduces the Bohr atom faithfully. 
Quantum mechanics in the Bohr atom is in fact the first approximation 
in the value v/c of the delay. 

The fact that the Planck constant is expressed in function of the 
parameters associated to the particular model system (see (|T2|)) could 
be a problem. If we consider another model system, usually another 
expression of the Planck constant would appear and another numerical 
value for the constant g would have to be chosen, and hence another 
expression for the delay. However, we think that in the expression 
of the delay, in fact, we will obtain the same constant g in all the 
problems as it happened in the development of the quantum mechanics 
(perhaps not exactly the same constant but a constant that will be 
expressible in terms of known constants and in function of the fine 
structure constant). Suppose, for a moment, that we are in 1913 and 
we made this work in that date. In order to be consistent with the 
Relativity theory, a new constant g needs to appear. This constant g 
seems to be another important constant of the nature. Suppose that we 
follow this line of research and we apply this theory to other systems. 
What we think that it would have happened is that in the resolution of 
these new systems, the constant g would also appear again and all the 
physicians would have thought that this g is a new important constant 
of the nature (in fact it is k 2 n/a). If we do not compare with the Bohr 
model, we find at the end that v/c = ir/(giri) and n/gi would be the 
relation between the velocity of the electron at the first energy level 
and the speed of the light (this relation is known in quantum mechanics 
as the fine structure constant a). 
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In the following section we will obtain a first approximation of the 
fine structure corrections of Sommerfeld j2U I2H] to the hydrogen atom. 




3. The fine structure of the hydrogen atom 

In old quantum mechanics, one of the more spectacular successes 
which helped to accept the Sommerfeld- Wilson-Ishiwara quantization 
was the study realized by Sommerfeld in 12*5] . These works treat 
with the hydrogenoid atoms giving an explanation of the fine structure 
of the hydrogen atom discovered by Michelson [T"~J E3] Sommerfeld 
applied Special Relativity theory assuming that the electron inside the 
atom is travelling near the speed of light. He had obtained the following 
correction to the energy levels 

(16) E^-^(Zac) 2 ^ 

where n is a positive integer number and n$ — 1, 2, 3, . . . , n. The agree- 
ment of this formula with the experimental data is fortuitous, because 
the correct formula is obtained using the Dirac equation, where the 
spin of the electron is taken into account. The final expression is the 
same that (|16j) changing by j + 1/2 where j is the total angular 
momentum of the electron and j can take the following values for the 
n level j = 1/2, 3/2, . . . , n - 1/2. 

In fact, the fine structure is a result of relativistic corrections to 
the Schrodinger equation, derived from the relativistic Dirac equation 
for an electron of mass m and charge e in an external electrical field 
— V$(r). Performing an expansion in v/c, the result for the Hamilton- 
ian H can be written as H = Hq + Hi, where 



2m Awe^r 

is the non-relativistic hydrogen atom, Z = 1, and Hi is treated as a 
perturbation to H , using perturbation theory. Hi consists of three 
terms: the kinetic energy correction, the Darwin term, and the spin- 
orbit coupling, Hi = H KE + H Darwin + ido- 
lising the delay theory introduced in this work we are going to ob- 
tain a first approach to the fine structure of the hydrogen atom, as 
straightforward application of the result of the previous section. 
The relativistic kinetic energy is given by Ex = /ic 2 — /ioc 2 , where 



/io is the rest mass and [i = /io/v 1 — v 2 /c 2 . Hence, the total energy is 
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given by 



E 



Hoc 



e 2 
K-, 

T 



^l-v 2 /c 2 

or (substituting the value of the Potential energy in terms of angular 
momentum) 



E 



Hoc 



y 7 ! - v 2 /c 2 



1 



K 2 hc 4 
L 2 



From the value of the angular momentum L = (Ke 2 k 2 )/(kac) of the 
previous section, we obtain 



E 



Hoc 



1 



yjl ~ V 2 /C 2 



1 



Hot 2 c 2 
k 2 



Substituting h — /Wyl ~~ v 2 /c 2 and developing in powers of v/c to 
fourth order we get 



E 



+ 



3^ 
4 c 4 



+ 



2 2 

Hoot c 
k 2 



f_ 3^ 
+ 2c 2 + 8 c 4 + " ' 



Finally, substituting the value of v — ac/k (obtained in the previous 
section) we arrive to the expression 

„2 



£^--Mac) 2 - 



1 + 



a 



1 - 



and, identifying n — \k\, it takes the form 



E 



1 



/i (ac) 



rr 



1 + 



a: 



n \n 



3 

4n 



If we compare this last expression with (|16p we see that we have ob- 
tained only a first approach to the Sommerfeld quantization because 
we only obtain the case = n. In this first approximation we have 
not introduced the intrinsic angular momentum of the electron (spin). 
Therefore further investigations are required to find the correct model 
of the delay equations with the intrinsic angular momentum of the 
electron taken into account. 



4. Concluding remarks 

The aim of my work is the understanding of which can be the expla- 
nation of the quantic behavior of nature. Physicians, working in quan- 
tum mechanics, believe that this explanation lies on the same theory. 
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However, there is no impediment to the existence of a deeper descrip- 
tion giving as an observable result the effects of quantum mechanics. 
In order to give an example, this is the same thing that happens with 
classical mechanics and statistical mechanics of n bodies. Under the 
laws of statistical mechanics, the laws of classical mechanics lie, that 
govern the movement of each one of the n bodies. This work does not 
pretend, in any case, to substitute quantum mechanics but to complete 
the knowledge that it gives to us. 

The whole foundation of quantum mechanics, starting from Heisen- 
berg and Schrodinger is somewhat ad hoc. This is not a great defi- 
ciency, since any theory must start with some assumptions, and the 
less is the number of such assumptions, the better is a theory. In this 
way, the model that we present in this work is better than the Bohr 
model because the unique assumption is that the electrodynamic inter- 
action has finite propagation speed. The introduction of some model, 
different from the existing one, should go accompanied with a possibil- 
ity of describing the new phenomena, yet to be discovered. We hope 
to give an answer to new phenomena in the context of this theory in 
future works. 

From the fact that the atomic Bohr model can be completely de- 
scribed by means of functional differential equations, we believe that 
in the future, observing the physical reality at a deeper level, we will 
be able to interpret the laws of probability of the quantum physics 
as the statistical results of values of certain variables perfectly deter- 
mined that at the moment are hidden for us... (we know that this idea 
has been extremely discussed) the history of science shows us that the 
current state of the knowledge is always provisional and that it should 
exist, beyond what is known, immense regions to be discovered. How- 
ever, functional differential equations are fundamentally different from 
ordinary differential equations, and their solutions can have qualita- 
tive features which are impossible for solutions of ordinary differential 
equations. The physical consequences of these differences are explained 
at length in |Tj|. We mean that the rich features that can appear in 
functional differential equations may explain quantum mechanics, from 
a deeper point of view. 
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